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ABSTRACT 

In order to precisely determine temperature and density of molecular gas in 
the Large Magellanic Cloud, we made observations of optically thin ^^C0( J = 3 — 
2) transition by using the ASTE 10m telescope toward 9 peaks where ^^CO(J = 
3 — 2) clumps were previously detected with the same telescope. The molecular 
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clumps include those in giant molecular cloud (GMC) Types I (with no signs of 
massive star formation), II (with HII regions only), and III (with HII regions 
and young star clusters). We detected ^^CO(J = 3 — 2) emission toward all 
the peaks and found that their intensities are 3-12 times lower than those of 
12C0(J = 3-2). We determined the intensity ratios of ^^00(7 = 3 - 2) to 
i3C0(J = 3-2), RfJ^^ and ^'^CO{J = 3 - 2) to ^'^CO{J =1-0), i?^^2/i-o' at 
45" resolution. These ratios were used for radiative transfer calculations in order 
to estimate temperature and density of the clumps. The parameters of these 
clumps range kinetic temperature Tkin = 15 - 200 K, and molecular hydrogen 
gas density n{B.2) = 8x10^ - 7x10^ cm~'^. We confirmed that the higher density 
clumps show higher kinetic temperature and that the lower density clumps lower 
kinetic temperature at a better accuracy than in the previous work. The kinetic 
temperature and density increase generally from a Type I GMC to a Type III 
GMC. We interpret that this difference reflects an evolutionary trend of star 
formation in molecular clumps. The RI^_2/i_q and kinetic temperature of the 
clumps are well correlated with Ha flux, suggesting that the heating of molecular 
gas n{B.2) = 10^ - 10^ cm~'^ can be explained by stellar FUV photons. 

Subject headings: galaxies: individual (LMC) — ISM: clouds — ISM: molecules 
— Magellanic Clouds — radio lines: ISM — submillimeter 



1. Introduction 



Star formation is of fundamental importance in understanding the evolution of galaxies. 
Stars are formed in dense clumps of giant molecular clouds (GMCs), and kinetic energy and 
heavy elements are ejected from stars back into the interstellar medium (ISM) through stellar 
winds and supernova explosions. This cycle enriches metal abundance in the ISM and drives 
the evolution of galaxies both physically and chemically. It is therefore important to better 
understand the evolution of GMCs, the principle sites of the star formation in galaxies. 

The Large Magellanic Cloud (LMC) is the most suitable galaxy to study star formation 
and natal GMCs because of its ideal location. The L MC offers a unique opportunity t o 
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achieve the highest resolution due to its proximity, 50kpc (iKeller fc Wood 
and nearly face-on position with an inclination angle of ~ 35° (Ivan der Marel &: Cioni 



19991) 



200111 



This also provides an advantage of less contamination along the line of sight. The LM C shows 



activ e on-going star formation and many massive young clust ers are being formed (IHodge 
1961 ). In the LMC, the metallicity is a factor of ~ 3 - 4 lower flDufourlll984l : iRoUeston et al.l 



2002) and the gas-to-dust ratio is a factor of ~ 4 higher ( Koornneei 1984 ) than values for the 
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solar neighborhood. The far-ultraviolet (F UV) radiation field is more intense in the LMC 
than in the Milky Way fjisrael et al.l Il986[ ) and the visual extinction is a few times lower 
than in the Milky Way. These infiuence the physical properties of GMCs, and may affect 
the initial condition of star and cluster formation. 

The first spatially resolved surveys of GMCs in the w hole LMC were made by the 
NANTEN 4m teles cope in the ^^CO{J = ] 



0) transition flFukui et al.lll999l 12001 



2008 



Mizuno et al.l 120011 ). These surve ys revealed t he dis tribution of the GMCs within a single 
galaxy at a ~ 40 pc resolution. iFukui et al.l ( l2008l ) derived the physical properties, such 
as size, line width, and virial mass, of the GMCs, and found that the ^^CO(J =1 — 0) 
luminosity and virial mass of the cloud s show a good correlation. Assuming that the clouds 
are in virial equilibrium, iFukui et al.l (|2008[ ) derived the X^o factor, conversion factor of 
the ^^CO(J =1—0) intensity to total molecular column density. The derived X^o factor 
was similar to that of other Local Group galaxies, such as Small Magellanic Cloud, M31, 
M33, ICIO, and Milky Way, and the CMC mass distribution, dN/dM, was also similar 
to that of M31, M 33, and ICIO, suggesting t hat GMCs in th e Local Group galaxies have 
similar properties (IFukui fc Kawamurall2010l : iBlitz et al.l 120071 ). Comparisons between the 
GMCs and signs of star formation such as HII regions and young stellar clusters were used 



to classify GMCs into three types in terms of star formation activities (IFukui et al.l Il999 



Kawamura et al.ll2009l ): Type I shows no signs of massive star formation. Type II associated 
with only small HII regions, and Type III associated with both HII regions and young stellar 



clusters, and these Types are interpreted as an evolutionary sequence (IKawamura et al. 
20091 ). Comparative studies of the ^^CO{J = 1 - 0) to HI ratio are also a key to understand 



the ev olution of GMCs and the CMC formation via conversion of HI into H2. IWong et al. 



(120091 ) made t wo dimensional, i. e., spatial, comparisons of the second NANTEN ^^CO(J = 
1 — 0) s urvev (jFukui et al.ll2008l ) with HI data set combined ATCA and Parkes Telescope 
surveys ( Kim et al. 12003 ). They found that significant HI column densities (> 10^^ cm~^) 
and peak brig htness tempe r ature s (> 20K) are necessary but not sufficient conditions for 
CO detection. iFukui et al.l (120091 ) compared the three dimensional data cubes, including a 
velocity axis in addition to the two spatial axes. They found that GMCs are associated with 
HI envelopes on scales of ~ 50 - lOOpc, and the HI envelopes may be in dynamical equilibrium 
or may be accreting onto GMCs to increase the mass via HI - H2 conversion. Subsequently 
a higher resolution (45" ~ 10 pc at 50kpc) survey of ^^CO(J =1 — 0) molecular cloud with 
the Mopra 22m telescope, the Magellanic Mopra Assessment (MAGMA), re vealed a more 
detailed distribution of the ^2qo( J = 1 - 0) emission in the individual GMCs (iHughes et al.l 



20101). iHughes et al.l (120 lOl ) presented that the physical properties of star-forming GMCs are 
very similar to the properties of GMCs without signs of massive star formation. 



The physical properties such as kinetic temperature and density of the molecular gas 
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in the LMC have been investigate d in the higher- J transitions {J 



J = 4 - 3, J = 7 - 6) of CO fe.g.. ISorai et al 



1999 



Israel et al 



2003 



Bolatto et all 2005: Kim et al 



2001 



ohansson et al 



2004 



Pineda et all2008l : lMizuno et all2010l ). Some authors (e.g. jBolatto et al 



Kim 



2006 



2 - 1, J = 3 - 2, 



1998 



iHeikkila et al. 



Nikolic et al 
2OO5I : 



2007 



Nikolic et al. 



20071 ) suggest two-component model especially for molecular clouds closely associated with 
HIT regi ons. Most of these stud ies targeted extraordinary HII regions so the sample is 
limited. iMinamidani et al.l ( l2008l ) have carried out ^^CO(J = 3 — 2) observations of 6 GMCs, 
including one Type I, two Type II, and three Type III GMCs, in the LMC with the ASTE 
10m telescope at a spatial resolution of 5 pc, and identified 32 molecular clumps. These data 
were combined with available ^^CO(J =1 — 0) and ^'^CO(J = 3 — 2) data and compared with 
LVG calculations for 13 clumps. The results show that these clumps range from cool (~ 10 
- 30K) to warm (> 30 - 200K) in kinetic temperature, and warm clumps range from less 
dense (~ 10'^ cm~'^) to dense (~ 10^'^ - 10^ cm~'^) in density, wherea s only lower lirnits in 
kinetic temperature were obtained in the warm clumps. Most recently, iMizuno et al.l (120101 ) 
have made ^^CO(J = 4-3) observations of the N159 region with the NANTEN2 4m sub- 
millimeter telescope. These data were used in LVG analysis combined with ^^CO(J = 1 — 0), 
(J = 2 — 1), (J = 3 — 2), and ( J = 7 — 6) as well as the isotope transitions of ^'^C0( J = 1—0), 
(J = 2 — 1), (J = 3 — 2), and (J = 4 — 3). The kinetic temperatures and densities were found 
to be ~ 70 - 80K and ^3x10^ cm'^ in N159W and N159E, and ~30K and - 1.6 xlO^ cm'^ 
in N159S, indicating that an analysis including higher- J transitions of both ^^CO and ^^CO 
can better constrain kinetic temperature and density. 

In the present study, we aim to determine temperature and density of molecular hydro- 
gen ,H2, gas in all three CMC types. Type I, II, and III, with higher accuracy by combining 
the ^^C0( J = 3 — 2) data newly obtained by using the ASTE telescope and ^^C0( J = 3 — 2) 
and ^^CO(J =1 — 0) data obtained with the ASTE and the SEST telescopes, respectively. 
All our data sets were convol ved to the s ame resolution of 45", corresponding to ~ lOpc at a 
distance of the LMC, 50kpc (lFeastlll999l ). and the large velocity gradient (LVG) calculations 
are performed to estimate the line intensities. 

Section 2 describes observations of ^'^CO(J = 3 — 2) transition. Section 3 shows the 
observational results and describes ancillary data sets. Section 4 shows data analysis, and in 
Section 5, we discuss the physical properties of clumps and evolutionary sequence of GMCs. 
In Section 6, we present the summary. 
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13 



CO(J = 3 — 2) observations 



2.1. Selection of molecular clumps 



The present targets are chosen from the 32 molecular clumps identified by the ^^C0( J = 
3 — 2) observation s of the six GMCs in th e LMC with the ASTE 10m telescope at a spatial 
resolution of 5pc (IMinamidani et al.l 120081 ) . The s ix GMCs we r e chos en from the NANTEN 
catalog of ^^00(7 =1-0) GMCs compiled by IPukui et all (120081 ). and they were three 
Type 111 GMCs, two Type 11 GMCs, and one Type 1 CMC. 

In the present study, we selected 9 molecular clumps in the Type 111 GMCs, LMC N 
J0538-6904 (the 30 Doradus region; 4 clumps), and LMC N J0540-7008 (the N159 region; 
3clumps), Type 11 CMC, LMC N J0532-7114 (the N206D region; 1 clump), and Type 1 
CMC, LMC N J0547-7041 (the GMC225 region; 1 clump). The observed clumps and their 
coordinates are listed in Table 1. Hereafter, the region names, which are parentheses above 
or column (4) in Table 1, and the numbers of clumps are used to identify clumps, as same 
as that in IMinamidani et al.l ( 120081 ). Among 9 molecular clumps, four clumps are in the 30 
Doradus region, 30 Doradus No. 1, No. 2, No. 3, and No. 4, three clumps are in the N159 
region, N159 No. 1, No. 2, and No. 4, one clump is in the N206D region, N206D No. 1, and 
one clump is in the GMC225 region, GMC225 No. 1. 



2.2. i'^CO( J = 3-2) observations 



The observations of ^'^CO(J = 3 — 2) transition at 330.587960GHz were made with the 
ASTE 10-m telescope at Pampa la Bola in Chile in Septem ber, 2006. In this term, a single 
cartridge type double-side-band (D SB) SIS receiver, S C345 (lKohndl2005l ). was installed and 
the XF-type digital autocorrelator (jSorai et al.ll2000l ) was operated in the wideband mode, 
providing a bandwidth of 512 MHz with 1024 channels. These correspond to a velocity 
coverage of 450 km s~^ and the resolution of 0.45 km s~^ at 330GHz. The half power beam 
width was 23" at 330GHz, and this corresponds to 5.6pc at 50kpc. We made 3x3 points 
mapping observations with position switching r nethod with a 20" grid sp acing toward the 
peaks of ^^CO{J = 3 — 2) clumps detected by IMinamidani et al.l (|2008[ ). Figure 1 shows 
the observed positions of the peaks of the ^^CO(J = 3 — 2) clumps named 30 Doradus No. 
1, No. 2, No. 3, No. 4, N159 No. 1, No. 2, No. 4, N206D No. 1, and GMC225 No. 1. 
The typical system temperature including the atmosphere was 260K in DSB. The pointing 
accuracy was measured to be better than 5" in peak to peak by periodically observing a CO 
point source R Dor in the ^^C0( J = 3 — 2) transition every 2 hours. The spectral intensities 
were calibrated by employing the standard room-temperature chopper-wheel method. We 
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observed Ori-KL and M17SW in the ^^CO(J = 3 — 2) transition to check the stabihty of 
the intensity cahbration, and the intensity variation during these observations was less than 
12%. 

The observed antenna temperature T^(OBS) varies with the side-band ratio (SBR) or 
the image rejection ratio (IRR): 



TX(OBS) 



TX(SSB) 
1 + 



IRR 



where T^(SSB) is the antenna temperature in the single side-band (SSB). We measured the 
antenna temperature, T^(OBS), and im age rejection ratio (IRR), simultene ously, with the 
newly installed 2SB receiver, CATS345 JSzawa et al.lboosi llnoue et al.lboosh . to the ASTE 
telescope in 2008, and derived the SSB antenna temperatures, T^(SSB), of IRC-I- 10216 
[R.A.(1950) = 09M5"'14.8% Dec.(1950) = +13°30'40"] and N159W [R.A. (B1950) = 05M0"03.7^ 
Dec. (B1950) = -69°47'00.0"]. From the observation of IRC+10216, the observed velocity in- 
tegrated intensity, /./.(T^(OBS)), was 56±4 K km s~^, and the image rejection ratio, IRR, 
was 6±2 (^8 dB) . This IRR was consistent with the results of the laboratory evaluation 
(llnoue et al.ll2008[ ). The /./.(T^(SSB)) was estimated to be 65±6 K km s~^. This value was 
compared with the data taken by CSO (jWang et al.lll994l ). and main beam efficiency at that 
moment was derived to be 0.70±0.06. From the observation of N159W, the observed antenna 
temperature, T^(OBS), was 1.83±0.04 K, and the image rejection ratio, IRR, was 4.7±0.4 
(~7 dB). This IRR was consistent with the results of the laboratory evaluation (Inoue et 
al. 2008). The SSB antenna temperature, T^(SSB), was estimated to be 2.22±0.06 K. We 
derive the main beam temperature of N159W to be 3.2±0.3 K. We scaled the observational 
data in 2006 to be consistent with this value. The achieved noise level ranges from 0.04 K to 
0.18 K in the main beam temperature scale. These observations were made remotely from 
NAOJ and NRO in Japan, by usi ng the network observation system N-C0SM0S3 developed 
by NAOJ jKamazaki et al.lboosh . 



Results &; Data sets 



3.1. i3co(J = 3-2) 



In Figure 2, ^^CO(J = 3 — 2) profile rnaps o f each observed peak are presented with 
12C0(J = 3-2) spectra JMinamidani et al.ll2008l ). The ^^CO{J = 3-2) intensities are 3 
- 12 times lower than those of ^^CO(J = 3 — 2). The peak velocities are consistent with 
those of i2co(J = 3 - 2), and the line widths are factor of 0.7 - 1.0 smaller than those 
of ^^CO(J = 3 — 2). Their line parameters at the peak positions of ^^CO(J = 3 — 2) are 
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summarized in Table 1. Detailed descriptions explaining the characteristics of each observed 
^^CO(J = 3 — 2) peak are given in Appendix. 



3.2. i2co(j = 3_2) 



We use the ^^CO(J = 3 — 2) data published by iMinamidani et all (120081 ) . which were 
taken by the ASTE telescope with an angular resolution of 22" . We measured the antenna 
temperature, T^(OBS), and im age rejection ratio (IRR), simulteneo usly, with the newly in- 



stalled 2SB receiver, CATS345 flEzawa et al.ll2008 



Inoue et al.ll2008l ). to the ASTE telescope 



in 2008, and derive the SSB antenna temperatures, TX(SSB), of IRC+10216 [R.A.(1950) 
= 09M5"'14.8^ Dec.(1950) = +13°30'40"] and N159W [R.A. (B1950) = 05M0°'03.7^ Dec. 
(B1950) = -69°47'00.0"]. From the observation of IRC-l-10216, the observed antenna temper- 
ature, r^(OBS), was 18.8±0.7 K, and the image rejection ratio, IRR, w as 22ib3 (^13 dB) . 



This IRR was consistent with the results of the laboratory evaluation (jinoue et al.ll2008l ). 
The SSB antenna temperature, T^fSS B), was estirnated to be 19.6±0.8 K. This value was 
compared with the data taken by CSO (jWang et al.lll994l ). and main beam efficiency at that 
moment was derived to be 0.60±0.02. From the observation of N159W, observed antenna 
temperature, T^(OBS), was 7.6±0.2 K, and image rejection ratio, IRR, was 9±2 (~10 dB). 
This IRR was consistent with the results of the laboratory evaluation (Inoue et al. 2008). 
The SSB antenna temperature, T^(SSB), was estimated to be 8.4±0.2 K. We derive the main 
beam temperature of N159W to be 13.9±0.7 K. We scaled the observational data taken in 
2004 ( IMinamidani et al.l 120081 1 to be consistent with this value. 



3.3. i3co(J=l-0) 



We use the CQ( J = 1 — 0) intensity of the 30 Doradus and N159 re gions published by 



ohansson et al.l (119981 ). and the N206D and GMC225 regions published by IMinamidani et al. 



(120081 ). These data are all taken by the SEST telescope, and the angular resolution was 45". 



4. Data analysis 

4.1. Derivation of line intensity ratios 

The spatial resolutions of the present CO data vary depending on the telescope and 
frequency. The angular resolutions of ^^CO(J = 3 — 2) and ^^CO(J = 3 — 2) data observed 
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by the ASTE telescope are 22" and 23", respectively, while that of ^^CO(J =1 — 0) data 
observed by the SEST telescope are 45". These correspond to ~ 5 pc and 10 pc, respectively, 
at 50kpc. We have convolved the ^^00(7 = 3-2) and ^^CO{J = 3-2) data into the 45" 
beam with a Gaussian smoothing function in order to derive physical properties of clumps 
whose sizes are around 10 pc. 

We fitted Gaussians to each of the ^^CO{J = 3-2), ^^CO{J = 3-2), and ^^CO(J = 
1 — 0) spectra with a single spectral peak. We derived the temperature peaks of the Gaussian 
curves, and the results are summarized in Table 2. The ratios of ^^CO(J = 3 — 2) to 
13C0(J = 3-2) (hereafter, i^J^^^) and ^^CO{J = 3 - 2) to ^^CO{J = 1-0) (hereafter, 
Rl^2/i-o) derived as the ratios of the peak values on a main beam temperature scale. 
The errors of RI^J2^ and -R35L2/i-o estimated to be 20 - 27% and 25 - 30%, respectively. 
A summary of these main beam temperature ratios and their accuracies are also presented 
in Table 2. These ratios will be compared with numerical calculations of radiative transfer 
in the LVG approximation to derive constraints on density and temperature . 



4.2. LVG analysis 

4-2.1. Calculations of a LVG model 



In o rder to estimate physica. 



properties of molecular gas, we have performed an LVG 
analysis flGoldreich fc Kwan Il974j) of the CO rotational transitions in the same way as de- 



scribed in lMinamidani et al.l (j2008l ). The LVG radiative transfer code simulates a spherically 
symmetric cloud of uniform density and temperature with a spherically symmetric velocity 
gradient prop ortional to the radius, and employs a Castor's escape probability formalism 
flCastorlll970l ). It solves the equations of statistical equilibrium for the fractional population 
of CO rotational levels at each density and temperature, because the level populations of 
CO are determined by both density and kinetic temperature of molecular hydrogen. It in- 
cludes the lowest 40 rotational levels of the ground vibrational level and use the Einstein A 
coefficient and H2 collisional ir npact rate coefficien ts obtained from the Leiden Atomic and 
Molecular Database fLAMDA: TSchoier et al1l2005[ ). 



We performed calculations of the fractional populations of the lowest 40 rotational levels 
of ^^CO and ^'^CO in the ground vibrational state over a kinetic temperature range of T^in 
= 5 - 200K and a density range of n(Il2) = 10 - 10^ cm~^. We assumed that the fractional 
abundance o f CO to H2 and the abundance ratio of ^^CO to ^^CO are 1.6x10"^ and 50, 
respectively jMizuno et al.lboiol : lilake et al.lll987h . 
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Figure 3 shows the general behavior of the loci of constant -R3_2/i_o "^•^'-'^ v^wiiDuant -itg.g 



and constant -Rs^^^^ 
is the ratio of transitions which have different 



in the density-temperature plane. The -Rg'ig/i-o 

critical densities, and traces excitation of molecular gas. The -^3%^^ is the ratio of transitions 
which have same critical densities and trace different column densities, and thus this ratio 
is a good tracer of column density and consequently volume density in wide temperature 
range. The combination of the two line ratios is nearly "orthogonal" , so the density and 



kinetic temperature are we l l cons trained. This is different from the case of iMinamidani et al. 
J2008h . IMinamidani et all J2008h used the combination of the ratios of ^"^COiJ = 3 - 2) to 
12C0( J = 1-0) (here after, i?3%/i_o) and ^^CO{J = 1 - 0) to ^'^CO{J =1-0) (here after, 
R^'^Jn^ ), which does not constr ain well the physical parameters for densities higher than 10^ 
cm~^ (IMinamidani et al.ll2008l ). 



4-2.2. Results of the LVG analysis 



We summarize the input parameters for the 8 clumps in columns (3) - (5) of Table 3. 
The higher transition data have been convolved into the 45" beam with a Gaussian kernel 
as described in section 4.1. The velocity grad ient, dv/dr, of each clump is derived from 
12C0(J = 3-2) data jMinamidani et allboosi l. 



Figure 4 shows the results of the LVG analysis for the 8 clumps. The horizontal axis 
is molecular hydrogen density, n{B.2), and the vertical axis is the gas kinetic temperature, 
Tkin- Solid lines represent -R3iL2/i_o ^^"^ dashed lines represent RI^J-^^ . The hatched areas 
indicate the overlap regions of these two ratios within the errors. In Table 3, we present the 
estimated densities and kinetic temperatures of the clumps. 



The derived density and kinetic temperature are well dete rmined compared to the pre- 
vious work based on the combination of Rl^2/i-o -^i-'o^^ ( Minamidani et al. 2008), es- 
pecially for the clumps in the warm Type III GMCs. T he results for the thre e clumps in the 
N159 region show a good agreement with the results of iMizuno et al.l (120 lOf ) who employed 
the calculations using high- J transitions of J = 4 — 3. 

Figure 5 summarizes the derived kinetic temperatures (T^in) and densities (n(H2)) of 
all 8 clumps. These range from 15 to 200 K in kinetic temperature and from 8x10^ to 
7x10'^ cm~^ in density. We found clumps in the 30 Doradus and N159 regions are warm and 
dense, except for N159 No. 4 (N159S), which is less dense and has an intermediate kinetic 
temperature. The clump in N206D, N206D No. 1, is dense and has an intermediate kinetic 
temperature. The clump in GMC225, GMC225 No. 1, is less dense and cold. The derived 
kinetic temperature and density denote the same tendency as indicated by Minamidani et 
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al. (2008) with improved accuracies. For example, the kinetic temperature and density 
of the N159 No. 1 (N159W) clump are determined with the acc uracies of 39% and 31% , 
respectively, in the present work, although, in the previous study (IMinamidani et al.ll2008l ). 
only lower limit of the kinetic temperature (30K) and the two orders of magnitude uncertainty 
of density (3x10^ - 8x10^ cm~^) were suggested for this clump. 



4-2.3. Discussion on the LVG results 



There are two distinctive aspects in the present work. One is the continuous coverage of 
molecular clump samples in the all types of GMCs, Type I, II, and III, which are interpreted 
as an evolutionary sequence of GMCs (Kawamura et al. 2009). The other is the combinations 
of line ratios used in the LVG analysis. The ^^GO{J = 3-2), ^"^CO^.] = 3-2), and 
^■^00(7 = 1 — 0) transitions are used for the analysis and the ^^CO(J =1 — 0) transition is 
not included in this work. This improves the accuracies of kinetic temperature and density 
of molecular clumps as described by Mizuno et al. (2010). 



W e also performed LVG calculations with the parameters used by IMinamidani et al. 
fl2008h . which are X{CO) = 3.0 xlO'^ and abundance ratio of ^^CO to ^^CO of 20, 25, and 



30. The kinetic temperatures and densities of clumps estimated from B^^_ll^ and R]^_2ii_q 



with these parameters were well determined. This indicates that the combinations of line 
intensity ratios used in this paper, RI^J2^ and Rl^2/i-o^ suitable to determine physical 
properties of molecular clumps in wide ranges. 



The N159 region is the most active site of high mass star formation in the Large Mag- 
ellanic Cloud. This has been observed in various molecular transitions by various telescopes 
located in the southern hemisphere, and many authors paid their effo rts to determine the 
physical properties and chemical coinposit ions via LVG analyses (e.g.. iMizuno et al.l 12010 



ul: 



Pineda et al.ll2008l : IMinamidani et al.ll2008l ) . In Table 4, some of the recent LVG results of the 
N159 region are summarized. The kinetic te niperature and density de rived in this work ar e 
consistent within errors with that derived by lMizuno et al.l (120 lOl ) and lPineda et al.l (120081 ) . 
which are determined with high accuracies based on data of higher transitions of J = 4 — 3 
and J = 7 — 6. 



4.3. Comparison of line intensity ratios and physical properties with Ha flux 

We use the H a data (Kim et al. 19991) t oward the present clouds using the method given 
in Appendix B in IMinamidani et al.l (120081 ). The typical background level of Ha flux is ~ 
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10~^^ ergs cm~^ at the 40" scale, which is a pixel scale of the Ha data. These data were 
regridded into the ^^CO(J = 3 — 2) data grids. The fluxes toward each ^^CO{J = 3 — 2) 
peak are listed in Tables 2 and 3. 

Figure 6 shows plots of R^^J-^^ and Rl^2/i-o functions of Ha flux. It is clear that 
-^3^2/1-0 ^^^^ correlated with the Ha flux with a correlation coefficient of 0.98, and there 
is no clear correlation between R^^ll^ and Ha flux. 

Figure 7 shows plots of molecular hydrogen densities, n(H2), and kinetic temperatures, 
Tkin, as functions of Ha flux. For dense clumps these are distributed over a wide range 
of Ha flux (10~^^ - 10"^*^ ergs s~^ cm~^), although lower density clumps are located only 
where Ha flux is weak (~ 10~^^ ergs s~^ cm~^). The Tkin of clumps are well correlated 
with the Ha flux with a correlation coefficient of 0.81. We note that the size scale here 
is lOpc and the good correla tion may not hold at smaller scales where local extinction of 



Ha becomes important (e.g.. iMizuno et al.ll2010[ ). These results suggest that far-ultraviolet 
(FUV) photons heat molecular gas whose density is 10^ - 10^ cm~^. This will be discussed 
in the following section. 



5. Discussions 

5.1. Evolutions of GMC and clump 

The results of our LVG analysis show that the clump kinetic temperatures, T^in, range 
from cool (~ 15K) to warm (~ 200K) and the densities, n(H2), from less dense (~ 8x10^ 
cm~^) to dense (~ 7x10^ cm~^). These large variations in the physical properties reflect the 
different characteristics of the GMCs according to their Types I, II, and HI, where Type I 
shows no signs of massive star formation. Type II associated with only small HII regions, and 
Type HI associated with both HII regi ons and young stellar c lusters, and these Types are 



interpreted as an evolutionary sequence ( iKawamura et al.ll2009l ). A clump in a Type I GMC 



GMC225 No. 1, has the lowest density (n(H2) ~ 0.99x10^ cm ^ ) and lowest temperature 
(Tkin ~ 25K). A clump in Type II GMC, N206D No. 1, is a dense (n(H2) ~ 3.2x10^ cm'^) 
clump and has intermediate temperature (Tkin ~ 42K). Clumps in Type HI GMCs, except 
for N159 No. 4 (N159S) clump, are dense (n(H2) ~ 4x10^ cm"^) and warm (Tkm > 50K) 
clumps. The physical properties are generally correlated with the star formation activity of 
GMCs. The Tkin and ri(H2) increase with the evolution of GMCs from Type I to HI. This 
trend becomes clearer with the present determination of Tkin and n(H2) than that of the 



previous study ( iMinamidani et al.ll2008l ). although the number of clump samples is smaller 
than that. 



-VI- 



The N159 No. 4 (NISQS) clump is a part of a Type III GMC, LMC N J0540-7008 
(iKawamura et al.ll2009[). This GMC is quite large and elongated from north to south (275pc 
X 53pc, P.A. = 87°; iFukui et al.ll2008l ). Young clusters and large HII regions are associated 
with the northern part of this GMC, and in the southern part, where N159 No. 4 (N159S) 
is located, young clusters and large HII regions are not associated with. The N159 No. 4 
(N159S) clump is less dense (n(H2) ~ 1.5x10^ cm~^) and has intermediate temperature 
(^kin ~ 39K). These p roperties are quite simil ar to that of the c l umps in Type I or II 
GMCs, as discussed by iMinamidani et al.l ( l2008l ) and iMizuno et al.l (120101 ). Comparison of 
the properties of molecular clumps with signs of star formation activities, such as young 
star clusters and Ha emission, at high spati al resolution will establi sh the classification of 
molecular clumps, whose sizes are ~ 7 pc (IMinamidani et al.l |2008[). Because the size of 
clusters in the LMC are distributed in a range from 0.1 - 10 pc (IHunter et al.ll2003[ ). which 
is smaller /similar size than the typical size of molecular clumps (~ 7 pc), it is important to 
compare molecular clumps with signs of star/clustar formation to understand the evolution 
of molecular clumps leading to the star/cluster formation. 



5.2. Heating of the molecular gas in the LMC 



As shown in the previous sections, the -Rglg/i-o correlated with Ha flux at a 10 pc 

scale, and the molecular gas kinetic temperatures (Tkin) are also well correlated with H„ flux. 
These present findings suggest that the heating of molecular gas whose densities are 10^ - 10"^ 
cm~^ may be dominated by far-ultraviolet (FUV) photons. The intense FUV field controls 
the physical and chemical processes in the ISM such as formation and destruction of molecules 
as well as ionization. These regions have b een modeled as photo-dissoc i ation regions (PDRs) 
or photon-domina ted regions (PDRs) (e.g., 
Rollig et al.l 120071 . and references their in). 



Tielens fc Hollenbach 1985: Kaufman et al. 1999 



The FUV flux [Gn) is estimated as 3500 in the 30 Doradus reg ion (IBolatto et al.lll999 



Poglitsch et al. 



19951: 



Werner et all Il978t Esrael fc Koornnee3 Il979f) an d 300 for the N159 



19961: 



region (iBolatto et al.l Il999l : llsrael et al. 

estimated to be (1 - 5) xlO"^ cm~^ for these regions. The PDR surface temper ature is 



Israel fc Koornneej|l979l ). The gas density is 

400 K 



for th e 30 Doradus region and ~ 200 K for the N159 region from Figure 1 of iKaufman et al. 
(119991 ). The PDR gas temperature is relatively constant from the cloud surface to a depth 
where either heating or cooling c hanges significantly. T he heating is generally dominated by 
the grain photoelectric heating (IKaufman et al.lll999l ). and then dust attenuation of FUV 
flux controls the thermal structure, where the typical size scale is 0.1 - 2 pc. These suggest 
that the effect of FUV heating of molecular gas seems to be local and direct phenomena. 
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The warm region can become larger under low-metallicity or high gas-to-dust ratio 
environments. These temperatures are basically consistent with the kinetic temperatures 
(^kin) of the warm clumps in the present sample if beam dilution by the present resolution, 
10 pc, is taken into account. 



5.3. Substructures inferred by the ^CO{J = 3 — 2) observations 



The ^^CO{J = 3-2) intensities are 3-12 times lower than those of ^^CO{J = 3-2) 
transition, and the line intensity ratios of ^^CO{J = 3 — 2) to ^^CO(J = 3 — 2), i?3^4^^. These 
vary not only from clump to clump but also inside of each clump. In the some clumps, clear 
two velocity components are detected in the ^^CO(J = 3 — 2) transition. These suggest some 
internal structures. 



Mizuno et all (120101 ) showed that the molecular distribution of the N159 No. 1 (N159W) 
clump in the ^^C0( J = 3 — 2) transition is similar to that of the i] Carinae northern cloud in 
the ^^CO(J = 1—0) emission line smoothed to 5 pc resolution. The original d ata of rj Carinae 
north ern cloud has a 2 pc resolution and several substructures are identified flYonekura et al. 
20051 ). This supports the existence of internal structures inside of ~ 5 pc scale molecular 
clumps, and indicates that these internal structures can be resolved in the CO transitions 
with high spatial resolution observations. 

Subsequent observations of t] Carinae northern cloud using a high density tracer such 
as H^^CO"'"( J = 1 — 0) at a high spatial resoluti on, resulted in the det ection of high density 
molecular cores whose sizes are less than 1 pc (lYonekura et al.ll2005l ). In the LMC, some 
initial interferometric observations of high density tracers, such as HCO"'" and NH3, were 
made with ATCA, using resolutions of Q" ~ ~ 19^^ correspondin g to ~ 1.4 - 5 pc at 50 kpc 
(e.g., IWong et al.ll2006[ lott eTahlboOsi boiol ). lott et all toO^ detected two peaks in the 



HCO"''(J =1 — 0) transition in the N159W region. Further systematic detailed observations 
using high density tracers at higher resolution, with ALMA for instance, should be important 
for probing the initial conditions of star/cluster formation. 



6. Summary 

We have made 3x3 points mapping observations in the ^^C0( J = 3 — 2) transition at 
330 GHz using the ASTE 10m telescope toward 9 peaks of ^^CO(J = 3 — 2) clumps, which 
cover all types of GMCs, Type I, II, and III, which are interpreted as an evolutionary 
sequences of GMCs. We have detected ^^CO(J = 3 — 2) emission from all peaks and 
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their intensities are 3-12 times lower than those of ^^CO(J = 3 — 2). From this, : We 
have derived the intensity ratios of ^^CO{J = 3-2) to ^^CO{J = 3-2), and 
^^CO(J = 3 — 2) to ^^CO(J =1—0), Rlt2/i-o^ resolution, and have compared these 

results with the LVG radiative transfer calculations in order to estimate kinetic temperatures, 
Tkin, and densities, n(H2), of the 8 clumps. The clumps show Tkin = 15 - 200 K and n(B.2) 
= 8x10^ - 7x10^ cm~^. The Ha flux toward these clumps is well correlated with the 
^^CO{J = 3 — 2)/^^C0(J =1 — 0) ratio, Rl^2/i-o^ ^^'^ with the kinetic temperatures, Tkm, 
of the clumps. Dense clumps are distributed in a wide range of Ha flux (10^^^ - 10~^° ergs 
s~^ cm~^), although lower density clumps are located only where Ha flux is weak (~ 10~^^ 
ergs s~^ cm~^). We found that clumps in the 30 Doradus and N159 regions (Type HI GMCs) 
are warm and dense, except for N159 No. 4 (N159S), which is less dense and has intermediate 
kinetic temperature. The clump in N206D (Type II GMC), N206D No. 1, is dense and has 
also intermediate kinetic temperature. The clump in GMC225 (Type I GMC), GMC225 No. 
1, is less dense and cold. We suggest that differences of these clump properties largely reflect 
an evolutionary sequence of GMCs and molecular clumps. The kinetic temperatures (Tkin) 
and densities (n(H2)) of molecular clumps increase generally with the evolution of GMCs 
and molecular clumps. The -Rgig/i-o kinetic temperatures (Tkin) are well correlated 
with Ha flux, suggesting that the heating of molecular gas whose densities (n(H2)) = 10'^ 
- 10^ cm~^ is dominated by FUV photons. The calculations of PDR models are consistent 
with this suggestion. 
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Promoted Research (20001003). T.M. was supported by JSPS Research Fellowships for 
Young Scientists. M.R. wishes to acknowledge support from FONDECYT(CHILE) grant 
N0IO8O335 and the Chilean Center for Astrophysics FONDAP No. 15010003. 

The ASTE project is managed by Nobeyama Radio Observatory (NRO), a branch of Na- 
tional Astronomical Observatory of Japan (NAOJ), in collaboration with Uversity of Chile, 
and Japanese institutes including University of Tokyo, Nagoya University, Osaka Prefecture 
University, Ibaraki University, and Hokkaido University. Observations with ASTE were car- 
ried out remotely from Japan by using NTTfs GEMnet2 and its partnet R&E (Research & 
Education) networks, which are based on AccessNova collaboration of University of Chile, 
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A. Detailed descriptions of each clump 

Detailed descriptions of cliaracteristics of eacti clump are presented in tlie following. 

30 Doradus No. 1 (30Dor-10) (Figure 2a) The measured noise level of ^^CO{J = 3-2) 
spectra in this region is 0.10 K rms at 0.45 km s~^ velocity resolution. The spatial extent 
and the velocity range of the ^^CO{J = 3 — 2) emission are well correlated to those of 
the ^^CO(J = 3 — 2) emission. The ratios of integrated intensities of ^^CO(J = 3 — 2) to 
^^CO(J = 3 — 2) (hereafter -^3%^//) are relatively small (~ 6) at the western position and 
relatively large (~ 19) at the northern position. 

30 Doradus No. 2 (30Dor-12) (Fig. 2b) The measured noise level of ^3C0( J = 3-2) 
spectra in this region is 0.09 K rms at 0.45 km s~^ velocity resolution. The spatial extent 
and the velocity range of ^^CO{J = 3 — 2) emission are well correlated to those of the 
^^C0( J = 3 — 2) emission. -Rs^^// at the northern position is ~ 3, the smallest value among 
this clump, and is larger than 15 at the southern position. 

30 Doradus No. 3 &: No. 4 (30Dor-06) (Fig. 2c) The southern part corresponds to 
No. 3 clump and the northern part corresponds to No. 4 clump. The measured noise level of 
^^CO{J = 3 — 2) spectra in this region is 0.08 K rms at 0.45 km s~^ velocity resolution. The 
spatial extent and the velocity range of the ^^CO{J = 3 — 2) emission are well correlated to 
those of the ^^CO(J = 3 — 2) emission. 

N159 No. 1 (N159W) (Fig. 2d) The measured noise level in this region is 0.18 K 
r.m.s. at 0.45 km s"^ velocity resolution. The spatial extent and the velocity range of the 
^■^CO(J = 3 — 2) emission are well correlated to those of the ^^CO(J = 3 — 2) emission. 
^3^-2^11 shows the smallest value of ~ 5 in this clump at the centeral position. At the south- 
west position, the -Rs^^// is also small (~ 6). At the north-west position, the ratio is large 
of ~ 15. 

N159 No. 2 (N159E) (Fig. 2e) The measured noise level in this region is 0.11 K 
r.m.s. at 0.45 km s~^ velocity resolution. The spatial extent and the velocity range of the 
^■^00(7 = 3 — 2) emission are well correlated to those of the ^^CO(J = 3 — 2) emission. 
^3^-2^11 shows the smallest value of ~ 5 in this clump at the western position and shows large 
value of ~ 15 at the north-east position. 



- 16 - 



N159 No. 4 (N159S) (Fig. 2f) The measured noise level in this region is 0.07 K r.m.s. 
at 0.45 km s~^ velocity resolution. The spatial extent and the velocity range of the ^^C0( J = 
3 — 2) emission are well correlated to those of the ^^CO(J = 3 — 2) emission. Intensities of 
the ^^CO(J = 3 — 2) emission are fairly weak compared to those of the ^^CO(J = 3 — 2) 
emission. R^l2n shows the smallest value of ~ 12 in this clump at the center position. 
Two velocity components are detected in the ^^C0( J = 3 — 2) transition, although only one 
component is detected in the ^^CO(J = 3 — 2) transition. 

N206D No. 1 (Fig. 2g) The measured noise level in this region is 0.09 K r.m.s. at 0.45 
km s~^ velocity resolution. The spatial extent and the velocity range of ^^CO(J = 3 — 2) 
emission are well correlated to those of ^2co(J = 3 - 2) emission. 

GMC225 No. 1 (Fig. 2h) The measured noise level in this region is 0.04 K r.m.s. at 0.45 
km s~^ velocity resolution. The spatial extent and the velocity range of ^^CO(J = 3 — 2) 
emission are well correlated to those of the ^^CO(J = 3 — 2) emission. Intensities of the 
^^CO{J = 3 — 2) emission are ~ 12 times weaker than those of the ^^CO(J = 3 — 2) 
emission. 
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Right Ascension (B1950) 

Fig. la.— Positions observed with ^^CO{J = 3 - 2) in the (a)30 Dor, (b)N159, (c)N206D, 
and (d)GMC225 regions. Observed positions are indicated with dots. Contours indicate 
integrated intensity of ^^C0( J = 3 — 2) and the grey scale indicates the Ha flux (Kim et al. 
1999). White crosses indicate positions of young cluster (<10 Myr; SWBO). 
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(a)30Dor No.1 
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Fig. 2a.— i3C0(J = 3-2) profile maps of the (a)30Dor No.l, (b)30Dor No.2, (c)30Dor 
N0.3&4, (d)N159 No.l, (e)N159 No.2, (f)N159 No.4, (g)N206D No.l, and (h)GMC225 No.l 
overlaid with ^^CO(J = 3 — 2) profiles. Black spectra are ^^CO(J = 3 — 2) and red are 
i2cO(J = 3-2). 
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(c) 30Dor No.3 & 4 
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Fig. 2c. 
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(d) N159 No.1 
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Fig. 2d. 
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(e) N159 No.2 
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Fig. 2e.— 
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Fig. 2f.— 
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(g) N206D No.1 
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(h)GMC225 No.1 
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Fig. 3. — Contour plots of LVG analysis for reference. Contours are (a) Rl^2/i-o' (^) -^3-2 5 
and (c) a + b. Here, X{CO)= 1.6 x 10~^, dv/dr = 1.0 km s~^pc~^, and the number ratio 
(abundance ratio) ^^cO/^^qo is 50. 
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Fig. 4. — Contour plots of LVG analysis for 8 clumps: (a) 30 Dor No.l, (b) 30 Dor No. 3, 
(c) 30 Dor No.4, (d) N159 No.l, (e) N159 No.2, (f) N159 No.4, (g) N206D No.l, and (h) 
GMC225 No.l. The vertical axis is kinetic temperature Tkin, and the horizontal axis is 
molecular hydrogen density n{B.2). Solid lines indicate Rl^_2/i~o' dashed lines indicate 
i?3%^^. Hatched areas are the regions in which these two ratios overlap within intensity 
calibration errors. 
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Fig. 5. — Plot of LVG results. The vertical axis is kinetic temperature, Ty^in, and the 
horizontal axis is molecular hydrogen density n(H2). 
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Fig. 6. — Relations between line intensity ratio and Ha flux, (a) Intensity ratio of ^^C0( J = 
3-2) to i3C0(J = 3-2), i^g^^l (b) Intensity ratio of i3qo(j = 3-2) to ^^CO{J =1-0), 
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Fig. 7. — Plot of physical properties as a function of Ha flux at ^^C0( J = 3 — 2) peak. The 
vertical axises are (a) molecular hydrogen density n(B.2), and (b) kinetic temperature, Tkin- 
The horizontal axis is Ha flux. 



Table 1. Observerd Clumps and Properties of ^^CO(J = 3 — 2) line. 





GMC 








Clump'^ 






Peak Properties of 


CO( J = 3-2) Line 






No.=' 


Name'' 


Type'' 






Position 


Ti-nh 


Vlsr 


AV 


Integrated Intensity 


rms 


Other ID'^ 








Region 


No. 


o(1950) 


5(1950) 


(K) 


(km s-i) 


(km s-i) 


(K km s ) 


(K) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


186 


LMC N J0538-6904 


III 


30 Dor 


1 


5 39 08.6 


-69 06 15 


0.99 


251.3 


5.5 


7.1 


0.10 


30Dor-10 










2 


5 38 54.6 


-69 08 00 


0.86 


247.4 


5.2 


5.2 


0.09 


30Dor-12 










3 


5 38 49.0 


-69 04 30 


0.54 


253.3 


3.8 


2.9 


0.08 


30Dor-06 










4 


5 38 49.0 


-69 03 30 


0.47 


249.0 


3.8 


2.2 


0.08 




197 


LMC N J0540-7008 


III 


N159 


1 


5 40 03.7 


-69 47 00 


3.2 


237.9 


5.8 


23 


0.18 


N159W 










2 


5 40 35.5 


-69 46 00 


1.7 


233.1 


3.7 


8.4 


0.11 


N159E 










4 


5 40 32.7 


-69 52 00 


0.60 


237.0 


3.5 


3.2 


0.07 


N159S 


156 


LMC N J0532-7114 


II 


N206D 


1 


5 32 58.4 


-71 15 20 


1.6 


225.4 


3.0 


5.7 


0.09 




225 


LMC N J0547-7041 


I 


GMC225 


1 


5 47 51.3 


-70 41 20 


0.19 


216.9 


2.2 


0.33 


0.04 





Note. — Units of right ascension are hours, minutes, seconds, and units of declination are degrees, arcminites, arcseconds. Col. (1): Running number of GMC used ioo 
Table 1 in Fukui et al. (2008). Col. (2); Name of GMC. Col. (3): Type of GMC. Col. (4): Region name used in this paper. Col. (5): Running number of ^'^CO{J = 3 - ^ 
clump in each region used in Table 2 in Minamidani et al. (2008). Cols. (6)-(7): Coordinates of the position of i2cO( J = 3 - 2) clump. Cols. (8)-(12): Observed propertied 
of the ^''CO(J = 3 — 2) spectra obtained at the peak positions of the ^^CO(J = 3 — 2) clumps. The peak main-beam temperature Tj^bi VlsRi and the FWHM line width 
AV are derived from a single Gaussian curve fitting and are given in cols. (8), (9), and (10), respectively. The ^^CO(J = 3 — 2) integrated inteisities and r.m.s. noise level 
at the peak positions of the '^^CO(J = 3 — 2) clumps are shown in cols. (11) and (12), respectively. Col. (13): Another identification based on ^^CO( J = 1 — 0) observations 
with SEST. 

"Fukui et al. (2008). 

''Kawamura et al. (2009). 

'^Minamidani et al. (2008). 

Johansson et al. (1998). 



Table 2. Line Intensities and Line Ratios. 



Region 
(1) 


No. 

(2) 


Convolved Tmb (K) 


^mb (K) 




Line Ratio 


Ha flux 
(xlO^-'-^ ergs cm^^) 
(9) 


i3cO{J = 3-2) 
(3) 


i2CO{J = 3-2) 
(4) 


i3cO{J = 1-0) 
(5) 


Ref. 

(6) 


^3-2 

(7) 


""3-2/1-0 

(8) 


30 Dor 


1 


0.47±0.08 


3.7±0.5 


0.13 


J98 


8.0±1.7 


3.6 ±0.9 


79 




2 


0.46±0.08 


2.5±0.4 




J98 


5.5±1.2 




59 




3 


0.30±0.06 


2.0±0.3 


0.13 


J98 


6.7±1.6 


2.3 ±0.6 


21 




4 


0.28±0.06 


2.3±0.4 


0.13 


J98 


8.2±2.2 


2.2 ±0.6 


10 


N159 


1 


1.5 ±0.2 


8.6±1.2 


0.80 


J98 


5.7±1.1 


1.9 ±0.5 


6.0 




2 


1.0 ±0.1 


6.1±0.8 


0.44 


J98 


6.1±1.2 


2.3 ±0.6 


9.3 




4 


0.36±0.06 


3.8±0.5 


0.72 


J98 


11 ±2 


0.50±0.13 


1.1 


N206D 


1 


0.86±0.12 


3.3±0.5 


0.85 


M08 


3.8±0.8 


1.0 ±0.2 


1.2 


GMC225 


1 


0.13±0.03 


1.5±0.2 


0.55 


M08 


12 ±3 


0.24±0.07 


0.94 



Note. — Col. (1): Region. Col. (2): Running number in each region. Cols. (3)-(4): The peak main-beam temperature, Tmb, of the ^''CO(J = 3 — 2) 
and ^■^CO(J = 3 — 2), respectively, derived by using a single Gaussian curve fitting for a spectrum obtaind by convolved spectra into the 45"beam with 
a Gaussian kernel. Cols. (5)-{6): The peak main-beam temperature, Tj^b! of the ^'^CO( J = 1 — 0) (col. (5)) and their references (col. (6)). Cols. (7)-(8): 
Ratios of the peak main-beam temperatures. Ratios of ^^CO{J = 3 — 2) to ^^CO{J = 3 — 2) are shown in col. (7), and ratios of ^^CO(J = 3 — 2) to 
'^^CO(J = 1 — 0) are shown in col. (8). Col. (9): Ha flux toward the peak positions. 



References. — (J98) Johansson et al. 1998; (M08) Minamidani et al. 2008. 



Table 3. LVG results. 



Region 


No. 


dv/dr 


pl2/13 


nl3 

3-2/1-0 


n(H2) 




Ha flux 






(km pc~^) 




(xlO^ cm-3) 


(K) 


(xlO^^^ ergs cm~^) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


30 Dor 


1 


0.9 


8.0±1.7 


3.6 ±0.9 


4.6 


134t^^ 


79 




3 


0.9 


6.7±1.6 


2.3 ±0.6 


4.5 tl;? 


78tf7 


21 




4 


0.5 


8.2±2.2 


2.2 ±0.6 


2.7 ±J-i 


13lt^^ 


10 


N159 


1 


0.9 


5.7±1.1 


1.9 ±0.5 


4.9 t\\ 




6.0 




2 


0.5 


6.1±1.2 


2.3 ±0.6 


o 7 +1.3 
•J.' -0.9 


lOOt^i 


9.3 




4 


0.4 


11 ±2 


0.50±0.13 


1 c; +0.3 
-0.2 


oq+25 


1.1 


N206D 


1 


0.3 


3.8±0.8 


1.0 ±0.2 


o 9 +0.9 
■^■^ -0.6 


42+15 

^^-13 


1.2 


GMC225 


1 


0.2 


12 ±3 


0.24±0.07 


".■^^-0.19 


2c: + 18 


0.94 



Note. — Col. (3): Velocity gradient of clumps (Minamidani et al. 2008). Cols. (6) and (7): Results of LVG 
analysis. 



Table 4. Summary of an excitation analysis in the N159 region. 







No.l (N159W) 


No.2 (N159E) 


No.4 (N159S) 


Paper 


Reference 


-^kin 


n(H2) 




^(H2) 


^kin 








(K) 


(10^ cm-3) 


(K) 


(10^ cm-3) 


(K) 


(10^ cm"3) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


This work 






A.9tli 


looli 


q 7+1.3 
' -0.9 


39tll 




MizunolO 


MIO 






oq+26 


q 1 +1.9 
'-'•-'--0.6 




-■-•"-0.3 


PinedaOS 


P08 


80 


10 - 100 










MinamidaniOS 


M08 


> 30 


3-800 


> 40 


1-300 


20-60 


1-6 



References. — 



(MIO) Mizuno et al. 2010; (P08) Pineda et al. 2008; (M08) Minamidani et al. 2008 



